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Abstract 

The increasing interest in compact astrophysical objects (neutron stars, binaries, galactic black 
holes) has stimulated the search for rigorous methods, which allow a systematic general relativistic 
description of such objects. This paper is meant to demonstrate the use of the inverse scattering 
method, which allows, in particular cases, the treatment of rotating body problems. The idea is to 
replace the investigation of the matter region of a rotating body by the formulation of boundary 
values along the surface of the body. In this way we construct solutions describing rotating black 
holes and disks of dust ( "galaxies" ) . Physical properties of the solutions and consequences of the 
approach are discussed. Among other things, the balance problem for two black holes can be 
tackled. 
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I. INTRODUCTION 



The systematic investigation of neutron stars and binaries consisting of pulsars and other 
compact objects and the increasing evidence for the existence of (rotating) black holes have 
stimulated theoretical and numerical studies on rapidly rotating bodies in General Relativ- 
ity. No doubt, realistic stellar models (e.g. neutron star models) require a careful physical 
analysis of their interior states and processes and, as a consequence, extensive numerical 
calculations. On the other hand, there is widespread interest for explicit solutions of the ro- 
tating body problem under simplifying assumptions. Such solutions could provide a deeper 
insight into physical phenomena connected with spinning matter configurations and, more- 
over, serve as test beds for the numerical investigations mentioned before. A good example 
is the Kerr solution, which has enriched our knowledge of rotating black holes in an in- 
estimable way. However, rigorous results for rotating bodies are relatively rare in General 
Relativity. Among other things, this is due to the mathematical difficulties with 'free bound- 
ary value problems', already known from Newton's gravitational theory, and to the specific 
complexity of the differential equations of Einstein's theory inside the body. Namely, the 
shape of the surface of a rotating self-gravitating fluid ball — the best model for astrophys- 
ical applications — is a 'compromise' between gravitational and centrifugal forces and not 
known a priori. (The surface is 'free', i.e. not flxed from the very beginning.) Though there 
are powerful (soliton-) techniques to generate (formal) stationary axisymmetric solutions of 
Einstein's vacuum equations, no algorithm to integrate the interior fleld equations is avail- 
able. Hence, at flrst glance, a boundary value description of rotating bodies seems to be 
questionable and inadequate. However, there are exceptional cases, in which the surface of 
the body has a known shape and the surface values provide enough information to construct 
the complete solution of the vacuum fleld equations. It is the intention of this paper to 
show that this is true for stationary black holes and disks of dust, which may be considered 
to be extremely flattened perfect fluid bodies. Moreover, it should become clear that our 
procedure, which is based on the inverse scattering method, opens an access to the not yet 
solved problem of the balance of two black holes and enables, in principle, the construction 
of black holes surrounded by dust rings ('AGN models'). 

Another interesting domain of application for the inverse scattering method is colliding 
gravitational waves. This theory is out of the scope of our paper. We refer to the article^ 
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FIG. 1: (a) A slice (j) = constant, t = constant, (b) 2-sheeted Riemann iC-surface with branch 
points (dots) and two cuts (sohd hnes)2i 

and references therein. 

The present work is mainly based on the papers^ and^ but it also contains substantial 
material not published before. 

II. THE BOUNDARY VALUE PROBLEM 

We consider a simply connected axisymmetric and stationary body and describe its ex- 
terior vacuum gravitational field in Weyl-Lewis-Papapetrou coordinates 

ds" = e-2^ [e^'^ (dp2 + dC') + pM02] - (dt + ad<pf (1) 

where the 'Newtonian' gravitational potential U and the 'gravitomagnetic' potential a are 
functions of p and C, alone. Fig. 1(a) shows the boundaries of the vacuum region: are 
the regular parts of the axis of symmetry (p = 0), S is the surface of the body and C stands 
for spatial infinity. Later on, we will integrate along the dashed line and pick up information 
from the boundary values of the gravitational fields at A^, B and C. The metric ^ allows 
an Abelian group of motions G2 with the generators (Killing vectors) 

= 51, < (stationarity) 

(2) 

?7* = 5^, Tj^rji > (axisymmetry) 
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where the Kronecker symbols 51 and 5^ indicate that ^* has only a t-component whereas r]'^ 
points in the azimuthal 0-direction (its trajectories are closed circles!). Obviously, 

e^^ = a = -e-2^r/,e (3) 

is a coordinate-free representation of the two relativistic gravitational fields U (generalization 
of the Newtonian gravitational potential) and a (gravitomagnetic potential). To get a unique 
definition of U and a, we prescribe their behaviour at infinity. Assuming that the space-time 
has to be fiat at large distances from the body and can be described by a Minkowskian line 
element in cylindrical coordinates, we are led to the boundary conditions 

C : U ^ 0, a ^ 0, k (4) 

Any linear transformation 

t' = t, (p' = <i)-ujt (5) 

introduces a frame of reference which rotates with a constant angular velocity uj with respect 
to that asymptotic Minkowski space. 

To describe stationarity and axixymmetry in that rotating system one would use the 
Killing vectors 

C = C'+u:^'\ V' = V'\ (6) 

instead of 

Regularity of the metric along means 

^± : a = 0, k = 0. (7) 

These conditions express the fact that is an axis of symmetry [a = 0) and ensure 
elementary flatness along the axis {k = 0). The behaviour of U and a at the surface B of 
the body depends on the physical nature of it. Rotating perfect fluids are characterized by 
a four-velocity field consisting of a linear combination of the two Killing vectors, 

= e"^ (f + n7]') , u'ui = -1 (8) 

where Q is the angular velocity of the body, and an invariant scalar pressure p, which is, for 
rigid rotation, 

Q = Qq, {Qq a constant) (9) 
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a function of V alone, 

p = piV), (10) 

as a consequence of the Euler equations. Along the surface of the body (if it exists) the 
pressure has to vanish, 

p(Vo) = 0, i. e. V must be a constant along B, 

B: e^^ = -(e + noV'm + ^o^) = e'''° (H) 

That is a further boundary condition. When identifying in (jSI) and VLq we introduce a frame 
of reference co- rotating with the body and may interprete V as the co-rotating 'Newtonian' 
potential, cf. and Interestingly, the event horizon 7i of a stationary (axisymmetric) 
black hole behaves like an 'ordinary' perfect fluid surface (jlip . Namely, one can show that 
a linear combination of the two Killing vectors, ^* + VL-aiff has a vanishing norm along 

H: e^"" = {C + ^im'm + ^m) = ^: (12) 

where VL^ is the angular velocity of the horizon. Hence we may include black holes in our 
scheme, see Fig. 1(a), for Vq — oo and H = B. It will turn out that (|T^ together with 
the correct positioning of the horizon Ti. in Weyl-Lewis-Papapetrou coordinates together 
with the asymptotic behaviour of the (invariant) potentials Q suffices for an explicit 
construction of the Kerr solution - thus providing a simple constructive uniqueness proof for 
stationary axisymmetric black holes. On the other hand, the condition (lllj) is not sufficient 
to calculate the gravitational vacuum field of rotating perfect fluid balls. However, in the 
disk of dust limit of such fluid conflgurations the fleld equations themselves will provide 
the missing boundary condition along the surface B of the disk, see^^. Starting with that 
completed set of boundary conditions we will be able to construct the global solution for 
the rigidly rotating disk of dust. 

III. THE FIELD EQUATIONS 

The vacuum Einstein equations for the metric coefficients k, U, a are equivalent to the 
Ernst equation 
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for the complex function 



where b replaces a via 



/(p,C) = e^^ + i& 



-4:11 U -iUl 

a^p = pe 0,^, a^^ = —pe b^^ 



(14) 



(15) 



and k can be calculated from 



= P 



Ul-Ul + ^e-^^bl-bl) 



(16) 



As a consequence of the Ernst equation (jlHj) . the integrability conditions a^p^ = a^^p and 
^,pC ~ k,Cp automatically satisfied such that the metric functions a and k may be cal- 
culated via line integration from the Ernst potential /. Thus, it is sufficient to discuss the 
Ernst equation alone. 



IV. THE LINEAR PROBLEM 

The existence of a Linear Problem (LP) for the Ernst equation^'^'^'i^ is the corner 
stone of our analysis since it provides a suitable instrument for tackling boundary value 
problems: the inverse scattering method (ISM). Here we will use a 'local' versioniA of the 
Linear Problem, 





0^ 




^0 




















^-4 





(17) 



where $(z, A) is a 2 x 2 matrix depending on the spectral parameter 




K-iz 

\ K + 1Z ^ ^ 

as well as on the complex coordinates z = p + i(,z = p — i(, whereas A, B and the complex 
conjugate quantities A, B are functions of z (or p, Q and do not depend on K. From the 
integrability condition and the formulae 
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it follows that a matrix polynomial in A has to vanish. This yields the set of first order 
differential equations 

= A{B -A)- ^(A + 5), = B{A - B) - ^(fi + A). (20) 

The system has the 'first integrals' 

A = B = (21) 

Resubstituting A and B in the equations (f^UI) one obtains the Ernst equation ()13|) . Thus, 
the Ernst equation is the integrability condition of the LP (fT7j) . Vice versa, if / is a solution 
to the Ernst equation, the matrix $ calculated from (fT7j) does not depend on the path of 
integration. The idea of the inverse scattering method (ISM) is to discuss for fixed but 
arbitrary values of holomorphic function of A (or K) and to calculate A, B and 

finally / afterwards. To obtain the desired information about the holomorphic structure in 
A, we will integrate the Linear System along the dashed line in Fig. 1(b) making use of the 
conditions (0]), ©, (fTTj) or (fT^ . In this way, we will solve the direct problem of the ISM and 
obtain ^{z,z,\) for z,z E A^^B^C. It turns out that the holomorphic structure remains 
unchanged by an extension of ^ off the axis of symmetry into the entire vacuum region such 
that one can construct functions $ with prescribed properties in A from which one obtains 
the desired solution /(z, z) everywhere in the vacuum region. This second step can be very 
technical and will, in general, lead to linear integral equations for In some circumstances, 
A may be replaced by K. For this purpose, it may be helpful to discuss the mapping ()18|) 
of the two-sheeted Riemann surface of K onto the A-plane for different values of p, C, (or 
equivalently z, z). Fig. 1(b) shows the position of the branch points = iz, Kb = —iz for 
the marked path A^CA^B of Fig. 1(a). It reflects the slice = constant, t = constant (Fig. 
^a)) and indicates, in particular, the position and shape of the body. Note that $ is not 
defined in the non-vacuum domain inside the circular contour around the origin. 

Consider now a Riemann surface with confluent branch points = = C ^ Here 
A degenerates and takes the values A = — 1 for K^s in the lower sheet, say, and A = +1 for 
K^s in the upper sheet {K ^ K^). 

We will now travel along the dashed line of Fig. 1(a) starting from and returning to any 
point p = 0, C ^ ■ (In Fig- 1(b) this corresponds to the bold faced points on the real 
axis.) Note that A = -1 for all iT's [K ^ Q in the lower and A = +1 for all K's {K ^ Q 



A Axis and Infinity 
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in the upper sheet of the Riemann X-surface belonging to axis values p = 0, ( E A"^ (the 
corresponding branch points cling to either side of the real axis in Fig. 1(b)). For p, C G C, 
the cut between the branch points (e.g., right solid line in Fig. 1(b)) points over the entire 
i^-surface and puts 'upper' K values into the lower sheet and 'lower' K values into the 
'upper' sheet. As a consequence, A will change from ±1 to =f1 between p = 0, ( = +oo 
and p = 0, ( = — oc^^. This 'exchange of sheets' is important for the solution of the linear 
problem: The initial value ^{pq,(oj ^) can (and must) be fixed only in one sheet of the 
i^-surface. The dependence on K in the other sheet follows by integration of the LP ()17|) 
along a suitable path^i. 

We will divide the integration of the LP (fTTjl along the closed dashed hne of Fig. [H^a) 
into two steps: 

(i) Integrating along A^CA^ 

This step can be performed without particular knowledge about the body and leads 
to a ^^geneml solution" for $ on the regular parts A^ of the symmetry axis. 

(ii) Integrating along B 

Here we confine ourselves to black holes and disks of dust. 

V. SOLUTION OF THE DIRECT PROBLEM 
A. Axis and Infinity 

Without loss of generality the matrix $ may be assumed to have the structure 



Note that both columns of $ are independent solutions of (|T7|). The particular form of (|^^ 
is equivalent to 




(22) 



together with 




(23) 




(24) 
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For K oo and A = — 1 the functions ifj, x may be normalized by 

^(p,C,-i) = x(p,C,-i) = i 

Finally, the solution to the Ernst equation can be read off at A = 1 {K oo) 

/(P, C) = x(p, C, 1), (7(M) = ^(P, C, 1)) • 



(25) 



(26) 



Remarkably enough, the Ernst equation retains its form in the frame of reference co-rotating 
with the body {u = Qq). This is a consequence of Q and (0) and implies the existence 
of a Linear Problem (|T7jl in the co-rotating system. In particular, the ^-matrices of both 
systems of reference are connected by the relation 



1 + iloa - ^]ope"^^ 








1 + Qoa + fiope"^^ J 



(27) 



+i{K + iz)noe 



-2U 



Henceforth, a prime marks 'co-rotating' quantities. We can now realise our programme and 
integrate the Linear Problem (fT7|) along the part A^CA^ of the dashed line in Fig. 1(a). 
Using (fT7|) along A^ and (pT| one finds for the axis values of $ 



A-" 




(28) 



A- : $ 



(29) 



/(C) 1 W 1 G{K) 

/(C) -1/ Vo 

where /(C) = /(p = 0, C) is the axis value of the Ernst potential and F{K),G{K) are 
integration 'constants' depending on K alone. The particular form of ()28|1 is due to the 
initial condition -i/^ = x = 1 for some po = 0, C = Co ^ -^^7 A = — 1 {K in the lower sheet), 
which fixes the second column of $ in (j^Hj) . cf. ()22|) . The first column corresponds to the 
upper (A = 1) sheet and represents a general integral with the two integration 'constants' 
F{K), G{K) which cannot be specified here. Along C, $ = ^{K) does not depend on p and 
C, since A and B vanish, cf. (PT|) . The 'exchange of sheets' along C, see Fig. Hfb), together 
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with leads to the particular form of $ on A~ . The representations (PHj). describe 
the behaviour of and x iii both sheets. Nevertheless, one may wish to consider the matrix 
$ as a whole as a unique function of A, which is therefore defined on both sheets of the 
i^-surface. From this point of view, the equations ()28j) . ()29|1 describe $ on one sheet only 
(say, on the upper sheet). Its values on the other (lower) sheet follow from ()24|1 . 
Combining ()28p .()29 p with ()27|) we obtain the axis values in the co-rotating system 



l + i{K- OQoe 



-1 -1 
1 1 



/(C) 1 FiK) 
/(C) -1 / I G{K) 1 



(30) 



A- 



l + i{K- OQoe 



-2U 



4 -1 
1 1 



X 



(31) 



/(C) 11 GiK) 

/(C) -1/ VO FiK) 

where 1 ist the 2x2 unit matrix. At the branch points = C of i^T-surfaces belonging to 
axis values p = 0, C G A^, ip X must be unique, i.e. 



A+iK^ = C) : * 



A-iK^ = C) : $ 





whence 



A 



A- 



+ . 



no 



no 



/(C) + /(()■ 



2/(0/(0 



^(0 



G(C) 



/(0-/(0 
/(O + /(O 



/(O - /(O 



(32) 



(33) 



(34) 



(35) 



/(0 + /(0 /(0 + /(0 

Thus, F{K) and G{K) consist in a unique way of analytic continuations of the real and 
imaginary parts of the axis values of the Ernst potential /(C)- Vice versa, /(C) follows from 
F{K),G{K) for K = (. Interestingly, the determinants of $ and can be expressed in 



B Surface 
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terms of ^f, K/' and F{K). From dHj) (Tr = (lndet$),^) , (EH) and (El -(EH), we 

have 

det $ = -2e2^F(ir), det = -2e2^F(ir) (36) 

where e^^ = 3?/ and e^^ = 3?/' (f/ = U{p, 0,V = V{p, ()). 

We may now interpret the result of (p?n|l - (p?H|) of the integration of the LP along A^CA~: 
On the regular parts of the symmetry axis, $ and can explicitly be expressed in terms 
of the axis values /(C) of the Ernst potential and its analytic continuations F{K),G{K). 
To calculate f{() one needs boundary values on B. Accordingly, the integration along B 
depends on the physical nature of the rotating body and can be performed in particular 
cases only. In the next section we will discuss black holes and rigidly rotating disks of dust. 

B. Surface 

1. One black hole 

We identify the surface B with the horizon Ti. In Weyl coordinates, the event horizon Ti 
of a single black hole covers the domaini^, 

H: p = 0, K^>C>K2. (37) 

(In Fig. ^a), the surface B degenerates to a 'straight line' connecting the regular parts 
A^,A'^ of the axis of symmetry.) Along H, e^^ has to vanish, see (EH), 

H: e'"" = {e + nov'm + ^om)=0 (l^o = ^^h)- (38) 

Because of 

cf. (H, this implies 

H: l + noa = 0. (40) 



B Surface 
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$ and can now be calculated along the horizon H. From dnZI),®, (001) and (|2Z|) 
we obtain 



/(C) l\ / U{K) V{K) 



JiO -Ij \WiK) X{K)^ 

/ -1 \ / f/(ir) 

= 2ifio(K-C) 

\ 1 / \ W{K) X{K) 

The Ernst equations have to hold at Ki and K2 too. Hence, $ and must be continuous 
in Ki and i^r2. Considering pH|) - (j!TT|) and (jlH), we are led to the conditions 

/i -1 \ / F 

/i + 2ifio(i^ - ^1) -1/ 1 

/i 



2ifio(^ - Ki) 




(42) 



/2 "M ( ^ ^ 

/2 + 2ifio(i^ - i^2) -1/ VO ^ 

/2 -1 

2ifio(^ - K2) 

where fi = f{( = Ki) and /2 = /(C = K2). Note that /i and /2 are imaginary, see (jSH). 
Eliminating the UVWX matrix, we obtain 




where 





(44) 



/F -G 

Ar= . (45) 

\G {1-G^)/f' 

Obviously, the elements of TV are regular everywhere in the complex i^-plane with the 
exception of the two simple poles at Ki and K2 {QKi = = QK2). The sum of the 
off-diagonal elements in ()45|1 must be zero. This requirement leads to the constraints 
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F{K) and G{K) take the form 

= Ani{K^ - Ki) ' = ^^m^-KiY ^^^^ 

Here we have chosen Ki = —K2, i.e., we have set the horizon in a symmetric position in 
the p, C-plane. Making use of (jH^ and and ehminating Qq by the second constraint 
equation we obtain the axis potential 

. f C(l + /i^) + (/i^-l + 2/i)i^i 

■ ^ c(i + /f) + (i-/? + 2/i)Kr ^^""^ 

It can be useful to introduce the multipole moments mass M and angular momentum J 
by an asymptotic expansion of /, 



and to replace /i , i^i in (gTj) , (gHl) and (gHl) : 

= K- + a--M^^ ^^^^ = K- + a--M- " ^^^^ 
To represent /(C), a simplifying parameterization is advisible, 



/i = itanv9/2, a = -Msnup, Ki = -K2 = VM^ - = Mcosip, ip = ^. (51) 
This yields 

. . _ (C-M) + iMsiny. 

^ • ^~ (C + M) + iMsin<^- ^ 

Finally, the second constraint equation becomes the well-known equation of state of 
black hole thermodynamics. 



2Ma,^^-jMl.,, (53) 

connecting the angular velocity of the horizon with mass and angular momentum. 



2. Two aligned black holes 

The same procedure can be used to tackle the balance problem for two black holes. 
The question is whether the spin-spin repulsion of two aligned stationary black holes can 
compensate their gravitational attraction. 
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Here we have two horizons Tii and 7^2 

n,: p = 0, K^>C>K2, ^2 : P = 0, > C > (54) 

separated by a piece of the regular symmetry axis 

^° : K2>C> Ks. (55) 

As a characteristic black hole property, the norm of the Killing vectors of the co-rotating 
frameworks has to vanish along the horizons, 

n,: ic, + nlr],){e + ^W) = o, {C, + nlv^){e + ^lvl = o, (56) 

where are the constant angular velocities of the respective horizons. 

Following the arguments for one black hole, we arrive ali^ 

where ^1 = ^2 = ^0' ^3 = ^4 = ^^"^ 

/ f. _i \ 

(58) 




whence 

FiK) - 



GiK) 



{K - Ki){K - K2){K - K^){K - K^) 

P2{K) 



(59) 



{K - K,){K - K2){K - Ks)iK - K,) ' 
where P4{K) and P2{K) are polynomials in K of the indicated orders. From (j59|) together 
with (jnH), (jHH)) we may read off the axis values of the Ernst potential. For the upper axis 
we obtain the structure 

m-^y (60) 

We need not use the representation of the explicit form of the second order polynomials 
?2(C)?Q2(C) cind of the constraints resulting from G = A/21 = —J^i2- Namely, from the fact 
that f{() is a quotient of polynomials of the same (even) order, it is clear that the desired 
two black hole solution can be generated by a Backlund transformation (in our case by a 
two-fold Backlund transformation) from the Minkowski space. (Note that the axis values of 
the Ernst potential determine solutions of the Ernst equation in a unique way.) 
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(a) 



infinity: 



disk: 



(b) 



12 



Kir 



FIG. 2: (a) Boundary value problem for the rotating disk of dust, (b) 2-sheeted X-surface of the 
rotating disk with branch points at — iz, iz^ 

The four constraints A/21 = —A/12 ensure that the constants Ki, fi = —f^ = 1, . . . ,4), 
Qq, Qq may be expressed by two position parameters and the masses and angular momenta 
of the two black holes. 

The Backlund generated solution belonging to (j^Oll known as the "double Kerr solution" 
was intensively discussed by several author o^'^i^^ . It turned out that there are necessarily 
struts between the "horizons". Since we have shown, by solving the boundary value problem, 
that Backlund generated solutions are the only candidates to describe aligned balanced black 
holes, we may now assert that black holes cannot be balanced at all. 



3. Rigidly rotating disks of dust 

Disks of dust can be considered to be extremely flattened spheroids (Fig. 2(a)) consist- 
ing of perfect fluid matter—. One can show that, for rigidly rotating dust, the boundary 
conditions and (jllj) have to be complemented by the condition b' = = on the 
disk (B). This condition follows from the Einstein equations as a transition condition from 
a divergence-free part of those equations via Gauss's theoren>^. Thus we have to take into 
consideration 



: regularity of /, B : f = e^^\ C : / ^ 1, (61) 
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see Fig. 2(a). On the disk, the Linear Problem of the co-rotating system of reference takes 
the form 

B : = ^ =^ I ° ^ 'M (62) 

where and /' are the 'co-rotating' $-matrix and the 'co-rotating' Ernst potential on the 
disk. This relation must be discussed under the 'boundary conditions' 



$'(p = 0, C = 0+, A) 1^ = $'(p = 0, C = 0+, A) 
cf. dSni) and 

$'(p = 0, C = 0-, A)|^ = $'(p = 0, C = 0-, A)|_^_, 

cf. daU- 

Again, this discussion allows the construction of F{K) and G{K) and, via (j34p and (|35|). 
the construction of the axis values f{() of the Ernst potential. We first take advantage of the 



symmetry of the problem which implies /(p, Q = /(p, — C) and connects the (^-derivatives of 
/' above {( = 0~^) and below {( = 0~) the disk 

TAc^o^ = -f'Ac^o- (63) 



A 



As a consequence, the LP ()62j) connects the matrix $ above the disk, 

A B B 

$= $'(p, C = 0+, K), with the matrix $ below the disk, $= $'(p, C = 0~, i^), 

A / l\ B 

B : ^= \ # H(ir), (64) 

-1 oj 

where the matrix ii{K) (the "integration constant") does not depend on p G i3. At the rim 
of the disk we have 

4(po,0,ir)=|(po,0,ir)=4. (65) 

_ r /O -l\ r 

Because of fl64j) . the rim matrix $ $ can be expressed in terms of 

V V 

A B 

$= $'(p, C = 0+,-?^), $= $'(p, C = ,K). Note that $ is considered to be a holomor- 
phic function of A and therefore a function living on the 2-sheeted Riemann i^-surface of 
Fig. ^b). Hence we have to discuss the rim matrix as a function of K on both sheets. 
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Any $ multiplied from the right by a matrix function of K is again a solution of the LP. 
The discussion of the rim matrix simplifies after the following redefinition 





Using (jnniKjSIl) we obtain 



where 



and 



e'^^'^AiS ^ on the upper sheet 
7^ = <( (67) 

-e^^'^S^^Ai on the lower sheet 



^ ^ G(i^) (G^ - 1)/f\ ^ ^ / /o/o - 4^]2ir2 ifeo + 2inoK ^^^^ 
-F{K) -G{K) / ' \ ^iVt^K -1 



/o = e2^° + i6o = /(C = 0+) . (69) 



/ -1 \ . /-I 

Note that A< = A/^ 

Obviously, Tr 7^ = Tr 7^ ^ = and ^ = 1 , whence 

Tr A<5-^ = Tr SM = 0. (70) 

This relation interlinks F{K) and G{K) and, because of (|34 |) . (|35| ) real and imaginary part 
of the axis values f{() of the Ernst potentia l^^i^^ . 

We next wish to determine F{K) and G{K) which in turn determine /(C)- To this end 
we consider $(p, C,A), for fixed coordinates p, C as a function of A. We have already used 
the initial conditions ip = x = for some p = po = 0, C = Co ^ •A.'^ prescribed in one sheet 
(A = —1) of the i^-plane. In principle the behaviour of $ in the other sheet and at all points 
in the p, C-plane can be calculated by integrating the LP along a suitable path. However, 
the coefficients A(p, C), B{p,() in the LP (jl7p are not explicitly known. Nevertheless, their 
regular behaviour outside the disk together with the boundary values on the disk, cf. (j^T|l . 
provides us with defining properties for One of them may be taken from Fig. ^b): Since 
the domain of the disk, < p < po, C = 0^ is a non-vacuum domain, where the LP fails, 
$ at the branch point pairs K = ip + 0^, — ip + 0^, < p < po cannot "pass" through 
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the contour T : — po < < po; i-G- ^ has a well-defined jump between opposite points 
along the contour F, see Fig. |21 A careful discussion would show that $(p, (, A), for fixed 
coordinate values p,( outside the disk (p, C ^ is a regular function in A outside F and 
jumps along F, i.e. ^» satisfies a (regular) Riemann-Hilbert problem. 

Consider now the jump where the signs indicate the two sides of F, cf. Fig. Efb). 

The LP tells us that does not depend on the coordinates and is therefore a function 

V of the contour alone, 

= Vi{K), K G Fi, (71) 

where u marks the upper and 1 the lower sheet. Since the jump contours Fu, Fi and the 
jump matrices V^, T>i are the same for all values of p, C (i-e. for all Riemann surfaces with 
different branch points), we may express and Vi in terms of the axis values of 

V,{K) =\'^ \ r ^ " I , KeV,. (72) 





A similar relation for Pi may be obtained via (j23)- As a consequence of (|7T| . ()72|1 the matrix 

(f q\ (f ^\ 

$ does not jump along Fu. Because of (12 7j) the same holds for and, 

finally, for TZ as defined in ()66|) . Consider now the Riemann K-surface of the disk rim p = po, 
( = 0. The cut between the branch points = ±ipo coincides with the contour Fu, F; 
which are on the two "bridges" connecting crosswise the upper with the lower sheet. Since 
TZ does not jump on Fu, we have, according to ^ {MS'^)- = -{S-^M)+. Though 7^ 
does not jump, F and G do jump, cf. (Tf^ . Note that F{K) and G{K) are unique functions 
of K. Hence, there is only one contour F : = 0, — po ^ < po where A4 (with the 
elements F{K), G{K)) does jump. Since $ is analytic outside Fu,Fi, the matrix M. must 
be analytic outside F. Thus we obtain F[K) and G{K) from the Riemann-Hilbert problem 

K : SM- = -M+S 

K ■ M{K) analytic in K, (73) 

S and Ai as in ()68|) . (Note that the elements of S, which are polynomials and the elements 
of which are rational functions in K do not jump along F.) There is no jump at the end 
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points of the contour K = ±ipo, A^(±ipo)- = ■M.{±ipo) + . As a consequence, one obtains 
TriS(±ipo) = 0, i.e., the parameter relation 

/o7o + Anlpl = 1. (74) 

It turns out that the Riemann-Hilbert problem (|7Sj) has a unique solution Ai{K) in the 
parameter region 

< /i = 2fi2e-2^0p2 < = 4.62966184 . . . (75) 

An important step on the way to this solution is the diagonalization of S. Finally, one 
obtains F{K), G{K) and the axis values of the Ernst potential f{() in terms of elliptic 
theta functions. We need not go this road. As we shall see in the next section, we can use 
the Riemann-Hilbert problem ()73j) to formulate a more general Riemann-Hilbert problem 
which will yield the complete disk of dust solution in terms of hyperelliptic theta functions. 

VI. ERNST POTENTIAL EVERYWHERE 
A. Kerr solution 

In the preceding section, we analyzed the axis values of the Ernst potential. We will 
now construct the complete solutions /(p, () of our boundary value problems from the 
information about the behaviour along the axis of symmetry gained by the discussion of 
the direct problem. 

There is, of course, no question that the discussion of the black hole case in IV B II will 
lead to the famous Kerr solution (in Weyl coordinates (P)). The point made here is that 
this solution describing the stationary rotating black hole can be derived from a boundary 
value problem. 

To achieve our goal it is useful to exploit the gauge freedom of multiplying $ from the 
right by an arbitrary matrix funktion of K. The transformation 

~ _ -a^ - fK + m ia \ 

preserves the properties (j22|) - (j25|) and enables the calculation of / via (j26|) . Because of p6|l . 
(I^mi. (pSj) and (fTH|). the determinant of $ becomes 

det$ = 7^^^|^(A^-A?)(A^-A^), 7 = 7(p, = - _ , (77) 



A Kerr solution 



20 



where 



This form of the determinant together with the axis values of $ tells us that $ must be a 
quadratic matrix polynomial of A, 

K -i- \7 

* = ^^(C + DA + EA2), (79) 
K 

where the 2x2 matrices C, D, E are functions of p, ( alone. It can be shown^*^ that $ with 
fl77|) satisfies the LP. (It is a Backlund transformation of the trivial solution / = 1.) 
According to (fTTj) . $(p, C, K) (^ = 1)2) must have a null eigenvector bi in the zeros Aj, 

#(p,C,A,)6, = (« = 1,2). (80) 

From the LP it follows that the elements of bi have to be constants. Hence, the quotient 

X(P,C,-A) C21-D2lA + ^2lA2' ^ ^ 

where the coefficients are elements of the matrices C, D, E must be a constant at A = Aj 
{i = 1,2). The values of the two constants (z = 1, 2) can be read off from the axis values of 
$ resulting from (jHH) together with (f^ . ^ and 

-icot-, = icot-. (82) 



X(-Ai) 2' x(-A2) 2' 

Note that x(~l) = 1 implies 

C21 - D21 + ^21 = 1. (83) 

These three conditions fix the coefficients C21, D21, E21 via a linear algebraic system. Finally, 
we obtain the Ernst potential everywhere from / = x(l)/x(— 1), 

_ rie^'^ + r2e'^-2Mcosv^ 

rie''^ + r2ei'^ + 2M cos 9? ' ^ ' 

where 

= (^.-C)' + p' (^ = 1,2) 



with Ki = —K2 and as in (|5T|1 . This is the Ernst potential / of the Kerr solution in Weyl- 
Papapetrou coordinates. By virtue of and this potential determines all metric 
coefficients in the line element (^. 
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B. Disk of dust solution 



In order to construct the ^-matrix for arbitrary values of p, ( and A, let us return to the 

^ , F ' 

Riemann-Hilbert problem ()7H|1 . As we have seen, the matrix $ does not jump 



G 1 



-1 

1 g' 



along Fu. Analogously, ^> does not jump along Fi. The images F^ of and T_x 



of Fi inherit these properties, which are essential to the following deductions. 

To formulate a Riemann-Hilbert problem in the A-plane, we define two matrices, 

1 



Q := e-2^o$ \{S + wl 

Fj \0 -1 




where 



Here we have made use of the parameter relation ()74|) . Since S and w are polynomials 
in K and therefore rational functions in A, the matrix Q has no jump at all. Taking the 
asymptotics of 5* and w into account, Q must take the following polynomial structure in A 

Q = iK + izf \ q^ = k\ + l\\ q2 = m + + p\\ q^ = q + r\^ + s\\ (87) 

yi3 -q2j 

where k, I, m, n,p] q, r, s are functions of p, ( alone. From the definitions (jHH|) . ()8(ij) and the 
condition (fTOj). we may derive 

QC = -CQ (88) 
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whereas the particular Riemann-Hilbert problem ()73|) has the continuation 

XeTx: (Q + e-^^^wl)C^ = -C+{Q + e-^^"^!) 

(Q-e-2^ou;l)/:_ = -/:+(Q-e-2^«wl) (89) 
A ^ Fa, r_A : analytic in A 

The following solution of the regular Riemann-Hilbert problem fl89j) is based on the 
diagonalization of Q. 

We consider a function \1/ defined by 



where 



+ e4Vb £22- Vl + «;2e-4yo£2,' ^ ^ 

^"1. (91) 

Note that \E' has no branch points at the zeroes Ki, K2, Ki = —K2 and K2 = —Ki of 

Kl = pf-^, Kl = pf-^ {^Ki < 0, > 0, /X as in ^) (92) 

(^ is unaffected by a change in the sign of yjw"^ + e'^^o). It is an odd function of A (vanishing 
at A = and at A = 00). Therefore, the function 

^ = ^/[A(fs: + i^)] = ^/^{K - iz){K + iz) (93) 

can be discussed as a unique function of K with the following properties: 

(i) Along r, because of (|H^. it jumps according to 

^- = ^+H In — ,. (94) 

(ii) Because of 

£^i(l + ^V^^')-£L = Q^„ (95) 

Q2I = - (i^ - i^a)(i^ - i^.), (96) 

the behaviour for K Ka/h is given by 

, ln(K-i^./fc) as (97) 

y/(ir,/,-iz)(ir,/, + iz)(^2^, + e4^o) 
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(The ambiguity of sign can be compensated for by the square root.) 

(iii) The behaviour for K oo, because of the definitions of Q and C, is given by 

as i^^oo. (98) 

These properties are reahzed by the following representation of ^: 



a 


- iz){K' + iz)y/w'^ + e^^o(^K' 
1 


-K) 




- iz){K' + iz)(w'2 + e4^o)(i^/ 
1 


I 




- iz){K' + iz){w''^ + e'^^o){K' - 


-K) 



-ipo 

-'I 

Kb 

- 2 / -w.., . ■ : . ... dir', (99) 

where Ka and have to be determined such that ^ = 0{K~^). The lower limits of 
integration in the last two integrals have been fixed to obtain the correct result in the 
Newtonian limit /i where Ka/Ki = 1 + 0{jj^) and Ki,/K2 = I + 0{jj^). (A systematic 
post-Newtonian expansion of the solution is given iniS.) Note that the last two terms in 
Eq. may also be interpreted as follows, 

Ka Kb Ki Kt Kt Kt 

2(y"+y")=2(y"{-}+y")=y"{i}+ y"{2}, (loo) 

Kl K2 Ka K2 Ka Ka 

showing that nothing special happens at Ki and K2. In this symbolic notation { — } indicates 
that the square root is meant to have the opposite sign with reference to the first term; {1} 
and {2} denote different paths in the complex i^-plane, which are chosen such that the 
closed integral 

Kt Ki K2 

j = j{2} = 2 j (101) 

Ka Ka Kl 

is performed around a contour enclosing the branch points Ki and K2 of Vw^ + e'^^o. In the 
subsequent formulae we normalize K and introduce 

X^!i, X^,,^^, X,,,-^. (102) 

Po Po Po 
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An asymptotic expansion of Eq. for X — oo {K oo) leads, according to (j^ . to 



In / = /i 



Xa Xi, 

x2 dx r x2 dx 
+ ■ w 



Xa Xi, 

dX /• dX 

J w 

X\ Xi 



w 

X2 

hdX 



hX^ dX 



Xa Xf, 

X dX /• X dX 

+ 



X-i X2 



hXdX 



(103) 



(104) 



where the lower integration limits Xi, X2 are given by 



2 1-/^ 



xl 



i + /i 



(3ftXi < 0, 3^X2 > 0), 



(105) 



H /i 

whereas the upper limits X^, X;, must be calculated from the integral equations (|104|) . Here 
we have introduced the abbreviations 

W = WiW2, Wi = v/(X - CIpoY + (p/po)^ W2 = ^l + ii\l + X''f (106) 



and 



In (v/l + /i2(l + X2)2 + /,(! + X2)) 



, . (107) 

7riv/l + p2(l+X2)2 

The third integral in ()103|1 as well as the integrals on the right-hand sides in ()104|1 have to be 
taken along the imaginary axis in the complex X-plane with h and and Wi fixed according 
to 'SiWi < (for p, ^ outside the disk) and 3?/i = . The task of calculating the upper limits 
Xa, Xh in (Uni from 

' hdX r hXdX 



u 



Wi 



(108) 



is known as Jacobi's inversion problem. Gopel^ and Rosenhain^l were able to express 
the hyperelliptic functions Xa{u,v) and Xfe(M, f) in terms of (hyperelliptic) theta functions. 
Later on it turned out that even the first two integrals in ()103|1 can be expressed by theta 
functions in u and f ! A detailed introduction into the related mathematical theory which 
was founded by Riemann and WeierstraB may be found inSSjiSiSi. The representation of the 
Ernst potential p03|) in terms of theta functions can be be taken from Stahl's book, see2^, 
page 311, Eq. (5). Here is the result: Defining a theta function {}{x,y;p,q,a) by 



00 00 



2mx+2ny+4:mna 



(109) 
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one can reformulate the expressions p03|) . ()lU4p to give 

^ ^(gpu + aiv - Ci, (3ou + (3iv - Ca; P, g, a) 
§{aQU + ait; + Ci, Pqu + (5iv + C2;p, g, a) 



with M and v as in ()108|) and 



^=1-^- (111) 



The normahzation parameters Oq, ai; /3o; Pi] 1o, li, the moduh p, q, a of the theta 
function and the quantities Ci, C2 are defined on the two sheets of the hyperelhptic Riemann 
surface related to 



W = ii^{X- Xi)(X - Xi)(X - X2)(X - X2)(X - i-z/p,){X + iz/po), (112) 

see Figure El There are two normalized Abelian differentials of the first kind 

dX XdX 
= ao— + ai—^ (113) 

^ dX ^ XdX 
dc.2 = /3o-^+A-^ (114) 

defined by 

dciJn = vri 5mn {m = 1,2; u = 1, 2) . (115) 



Eq. (jllSp consists of four linear algebraic equations and yields the four parameters ao, ai, 

Po, Pi in terms of integrals extending over the closed (deformable) curves ai, 02- It can be 

shown that there is one normalized Abelian differential of the third kind 

dX XdX X^dX 
d^ = 7o-^ + 7i-^ + /^-^ (116) 

with vanishing a— periods, 

^dw = (j = l,2). (117) 



This equation defines 70, 71 (again via a linear algebraic system). The Riemann matrix 

(In p 2a \ 
(t= 1,2; J = 1,2) (118) 
2a \nq J 

(with negative definite real part) is given by 

Bij = i dujj (119) 
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FIG. 3: Riemann surface with cuts between the branch points Xi and Xi, X2 and X2, — i-^/po 
iz/pQ. Also shown are the four periods and bi {i = 1,2). (Continuous/dashed hues belong to 
the upper/lower sheet defined by — > ±/nX^ as X ^ 00.)^ 

and defines the moduli p, g, a of the theta function ()109|) . Finally, the quantities Ci, C2 
can be calculated by 

CXD + 

Q = - y dcu, (« = 1,2), (120) 

-i^/po 

where + denotes the upper sheet. Obviously, all the quantities entering the theta functions 
and the exponential function in (jllOj) can be expressed in terms of well-defined integrals 
and depend on the three parameters p/po, C/Po, P- The corresponding "tables" for a^, (3i, 
7j, Cj, -Bjj, M, f , w can easily be calculated by numerical integrations. Fortunately, theta 
series like (jl(J9p converge rapidly. For < p < po? the solution (jllUj) is analytic everywhere 
outside the disk — even at the rings — iz/po = Xi, X2. The complete metric, calculated 
according to ((T)) and (fT ^ - (fTBj) is given in the Appendix. 

In the framework of the completely integrable evolution equations, the solution ()11U|I 
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FIG. 4: Relation between 2QqM and M"^ / J for the classical Maclaurin disk (dashed line), the 
general-relativistic dust disk and the Kerr black-holei^. 

may be interpreted as a 'Backlund-like' transformation of well-defined 'seed' solutions 
u,v,w satisfying axisymmetric Laplace equations. The transformation 'parameters' 
ao, /3o; tti, A; 70, 7i; Q', a; C15 C*2 depend on the 6 branch points of the 2-sheeted Riemann 
i^-surface associated with the function W = W{X), cf. ()106|) . and do not depend on u, v, w. 
All in all, / is a function of the 2 parameters po and fi and the 2 cylindrical coordinates p 
and (. For p ^ 1 we obtain the Maclaurin disk as the Newtonian limit. 

VII. PHYSICAL DISCUSSION 

Since the Kerr black hole is also a 2 parameter solution it might be interesting to compare 
the behaviour of both solutions in dependence on common parameters, say, on mass M and 
angular momentum J. It must be possible to express the area of the horizon and the disk, 
the radius po of the disk or other physical quantities in terms of M and J. A very illustrative 
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(a) 1-M/Mo (b) 1-M/Mo 

FIG. 5: (a) ^qMq and (b) J/Mq as functions of the relative binding energy (Mq — M)/Mq for the 
disk of dusti^ 

relation is the angular velocity VLq as a function of M and J, since VLq is defined in both cases. 
For black holes, we have derived the explicit expression Surprisingly, the corresponding 
disk of dust relation has the same scaling behaviour, i.e., MVLq is a function of M"^ / J alone. 
Fig. m shows this dependence for both solutiong^'^. For M'^ / J — ^ 1 (corresponding to 
/i — > /io), where the disk solution becomes identical with the extreme Kerr solution outside 
the horizon (p^ + > 0); there is a "phase transition" between the disk and the black 
hole. Note that for non- vanishing f2o, Po as /i ^ Po- A detailed analysis of the disk 
solution for /i ^ po? including the discussion of a different, non- asymptotically flat limit of 
space-time, which is obtained for finite p/po and C/po (p^ + = !), can be found in^^. 

We remark that (jllUj) solves the Bardeen- Wagoner problen>2^ explicitly. All metric coef- 
ficients in (1) are analytic in p, Q outside the disk and continuous through the disk. From 
a physical point of view we have an extremely flattened rigidly rotating body and, likewise, 
a rotating continuous distribution of mass points interacting via gravitational forces alone 
('galaxy' model). Fig. El illustrates the 'parametric' collapse of a disk with the total mass- 
energy M, the baryonic mass Mq, the angular velocity VLq and the angular momentum J 
towards the black hole limit (1 — M/Mq = 0.3732835 . . .). Imagine a disk consisting of a 
fixed number of baryons (fixed Mq): Occupying states with decreasing energy M, it would 
shrink thereby shedding angular momentum but increasing its angular velocity. The above 
mentioned limit of the relative binding energy 1 - M/Mq corresponds to the extreme black 
hole limit. Additional physical effects (ergozones, dragging effects, surface mass density . . . ) 
as well as further parameter relations have been discussed in^ and^i. 
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The methods outhned in this paper could be used to construct self-gravitating disks 
around a central black hole. 

APPENDIX 

The metric functions e^^, a, e^^ calcutated from the Ernst potential via (fT^ - ()lfi|l 

are given as follows: 



'i9(0)^9*(0)^9(a + c)i9*(a + c 

(1 

+ 



^ {l + noa)e^^ _ '(?(0)r(0)^9(a + 2c)r(a 



ilop ??(c)?9*(0)^9(a + c)?9*(a + c)^ 

2fc(p,C) _ ^(P' C) 



e 



with 



^9(a)^9*(a 



a) /, 9Mn^9(x)r(x) A 



where 



1 1 



2ko = ^ I fdXdX'^^'^'^ h{X)h{X'){X - X,){X - X,){X' + X,){X' + X,) 



A J J XX' (X-X') 



X + iz/po yX' + iz/po' 

?9(x) = 7?(x; p, g, a) = Xs; p, g, a), 

ivr Itt 

■(9*(x) = i3{xi + —,X2 + —■,P,q,a), 

a= {ai,a2) = {aoU + aiV,(3oU + Piv), 0=(0,0), c=(Ci,C2). 
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